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I . INTRODUCTION 



The program evolved during this study models a submerged 
circular ring using trigonometric series. It was developed 
for use in conjunction with a finite element fluid model 
based on a displacement potential formulation. The purpose 
of the combined models is to predict the effects of 
cavitation on underwater shock loading of the structure. 

More information on the fluid model can be found in 
References 1 and 2. The purpose of this paper is to 
describe the development of the structural model and to 
present some of the results obtained from its use in 
combination with the fluid model. A listing of the FORTRAN IV 
program implementing the structural model is given in the 
Appendix. 
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II. GOVERNING EQUATIONS 



A. DIFFERENTIAL EQUATIONS FOR THE DEFLECTION 
OF A THIN CIRCULAR RING 

In Figure 1 the solid line represents the ring after 
deformation, and the dotted line the ring before deformation. 
For small deflections the curvature of the element m^n^ can 
be taken as 



_1_ _ d9 -I- Ad6 
R^ ds + Ads 



where R, w, 0 and s are defined as shown in Figure 1 and 
w is taken positive inward. Making use of the relations 



Ade = de 1 - de 





ds 



dw 

ds 




ds 



Ads = ds^ - ds = (r-w)d0 - rd9 = -wd0 



-w 



ds 

R 



and substituting into the equation above yields 



_1_ 

R, 



.2 

d9 + — j ds 
ds 

dsd-^) 



d9 (1 + |) 
dsU-|) (1 + g) 



5%dsa + |, 

ds 

ds(l-g) (1 + |) 



Neglecting the higher order terms, the above expression 
becomes 
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/ 




Figure 1. Geometry of Ring Deformation 
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-L = de ( w 

R 1 ds 1 R ; 



d 2 w 

ds 2 



I (1 + 5) 



,2 

d w 

, 2 
ds 



or alternately 



1 _ _ 1 
R 1 ~ R 





(w + 



A 2 
d w, 



d6‘ 



( 1 ) 



For a ring where the thickness is small compared to the 
radius and elastic behavior is assumed, it can be shown 
that the approximate relationship between deflection and 
loading is [Reference 7] 




where M is the bending moment about the centroidal axis 
and D is the flexural rigidity of the ring. A positive 
bending moment produces compression in the outside fibers 
of the ring. Combining equations 1 and 2 yields the 
differential equations for the deflection of the ring given 
below. 




« 2 

. 3 w 



36 ‘ 



+ w) 



M 

D 



(3) 
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B. GOVERNING EQUATIONS OF MOTION 

The governing equations of motion used in the program 
were arrived at by the application of Hamilton's principle. 

In order to apply Hamilton's principle, it is first necessary 
to derive expressions for the strain energy and kinetic 
energy of the ring as well as the work done by the external 
loads. In these derivations the ring was taken to be 
homogeneous, elastic, and of unit width. The pressures on 
the ring and its deflection were represented by the 
pressures at and the deflection of a set of nodal points 
equally spaced along the circumference of the ring. 

. The shock front is assumed to approach the ring normal 
to the 9=0 plane as shown in Figure 2, where 9 is taken 
to be positive counterclockwise. 

K 



Ur 




Figure 2. Ring, Shock Front Orientation 
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Since the ring is symmetric about the 9=0 plane, 
results are given for 9 from 0° to 180°. The response of 
the ring is represented by the response of a set of N+l 
nodal points equally spaced along the ring. The radial 
displacement is approximated by a trigonometric cosine 
series 

N 

w = l a cos n9 . (4) 

n=0 n 

The tangential deflection is represented by the trigonometric 
sine series 



N 

v = / b sin n9 , (5) 

L 1 n 
n=l 

where v is taken to be positive in the negative 9 direction. 
Similarly the normal pressure applied to the ring is 
represented by 



N 

p = l c cos n9 . (6) 

n=0 n 

1 . Strain Energy 

The strain energy is comprised of two components 
[Ref. 4] . One component is the strain energy due to bending 
and the other to strain in the 9 direction (eg). The strain 
energy stored in the ring due to bending is calculated as 
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the work done by the internal moment acting over the 
entire ring, and for a ring of unit width is 



U. 



B 



/ M RKd6 

0 



where < is the curvature of the ring. Combining this 
equation with equation 3 and using the relation k = M/D 
yields 



The strain energy due to the strain in the 9 direction is 
obtained as the distance traveled by the average normal 
force in the 9 direction. 

The circumferential strain is composed of two 
components the first of which is the result of the radial 
displacement of the ring ( ) and the second results from 

I\ 

tangential displacement and is (^- . The strain energy 

K do 

from the normal force (N) is 



U. 



B 




(7) 



U 



9 



1 

2 



2ir 

/ Ne Q R d9 
0 



1 

2 



2tt 

/ Ae g E £ q R d9 

0 



or 



( 8 ) 



U 



9 



J / RAE£ 
^ 0 



9 



2 



d9 = T / (fj + w) 2 Rd9 

2 0 R 2 39 
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where A is the cross sectional area of the ring and 

E is the modulus of elasticity of the ring material. 



Combining equations 7 and 8, the expression for the 
total strain energy is 



U n 



, 2tt ,2 

1 f D , 8 w 

o ; "3 



2ir 



80' 



■. , Q . 1 ( AE , 8v ,2 JQ 

+ w) d8 + J + w) d0 

0 (9) 



Substituting into this equation the series expressions for 
radial and tangential displacements, and assuming constant 
geometrical and material properties over the ring yields 
after integration 



U„ 



ttD 

2R 3 



N 

I 

n=0 



2.2 . ,2 , AEtt 

a n (n - 1} + TF 



N 

l 

n=0 



(nb 



n 



+ V' 1 



An expression for the change in strain energy SU 
corresponding to small changes in tangential and radial 
displacements (6h^and6a n ) can now be found 

m N 

ttD r r , 2 . . 2 » , , AEtt , v , 2, , , 

5U = — r [ 2, (n - 1) a 6a + — — [ ) (n b^ + na ) Sb 

3 n=0 n n R n=0 n n n 

+ (nb n + a n ) 6a^) ] (10) 



★ 

In this equation and several that follow which contain 
a summation expression starting with n = 0, a factor of 2 
associated with the n = 0 term is omitted for the sake of 
brevity. The term is accounted for in the solution process 
used with the model. 
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2. Work 



The work associated with a small change in radial 
displacement is given by 



2tt 

6 W = / p 6w R d0 

0 



Substituting the relations 



N 

p = 7 c cos m9 

n m 

m=0 

and 

N 

6 w = ? 6a cos n8 

L n n 



into this expression yields, after integration 



N 

6W = y c 6 a ttR 
L n n n 



3 . Kinetic Energy 

The expression for the kinetic energy of the ring 
can be arrived at by considering an infinitesimal section 
of the ring that has been set in motion. The kinetic energy 
of the element can be represented as the sum of a contribu- 
tion due to translation of the mass center and a contribution 
due to rotation. This leads to 
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dT = j p (ARd0 ) (w 2 + v 2 ) + j pIRd6(^) 2 



where I is the area moment of inertia about the centroidal 
axis . 



Using the series relations 



N 

w = y a cos n6 
n=0 n 



N 

w = 7 a cos n6 

L n n 

n=0 



N 

v = J b sin n0 
L , n 
n=i 



N 



v = 



y b sin n6 
L n 



n=l 



yields 



1 N • 2 N 2 

dT = i p (ARd6 ) [ ( l a_ cos n6 ) + ( £ b sin n6 ) ] 

Z n=0 n=l n 

1 , N - 

+ -j pIRd9 (- =• E a n sin n6) 
n=l 



Integration over the ring yields 



N 



T = ^ l [(AR + n 2 |)i n 2 + ARb n 2 ] 



n=0 



The resultant change in kinetic energy due to a 
small change in velocity components can then be found as 



N 



6T = ttp E [ (AR + n 2 |)a n 6a n + ARb ob n ] 
n=0 



( 12 ) 
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4 . Application of Hamilton’s Principle To Find 
Coupled Equations of Motion 

From Hamilton's principle [Ref. 5] it is known 

that 

/ (ST - 6U + 5W) dt = 0 (13) 



Combining equations 10, 11, 12 and 13 yields, after 
carrying out an integration by parts on the first term. 



N 



l tt p ( AR + n 2 h [a 6a 
^ R n n 

n=0 



f a 6a dt] 
J n n 

r l 



N 

+ y it pAR [b 6b 

n=0 n n 



2 .. 

/ b 6b dt] 
t n n 



~ / 2 ( I | [ I (n 2 -1) ^ a Sal + 



N 



2 , \ 2 



AEtt 



ti R n=0 



n n 



N 2 

[ 7 (n z b + na ) 6b 

L n n n n 

n=0 



N 

+ (nb + a ) 6a ] - > c 5a ttR} dt = 0 

n n n „£ n n n 
n— u 



Since the 6a^ and 6b n can be chosen arbitrarily/ this 



becomes 

( i + <|> 



where c 



2 2 . - 
n )a n + 



b + (§) 2 
n R 



= (E/p) 



1/2 



<f ) 2 [1 + (|) 2 (n 2 -l) 2 ]a n + (|) 2 = 

na n + <l> 2 n \ * 0 

is the speed of sound in the ring 



c 

n 

PA 



(14) 



and 
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where r is the centroidal radius of gyration of the ring 

2 

section and the relation I = Ar has been used. 

5 . Separated Equations of Motion 

It is common practice to separate ring deformations 
into flexural modes and extensional modes. The initial 
version of this program utilized such a resolution. The 
flexural mode is characterized by the requirement that 

3v 

£ 0 = - (w + |qO/R = 0 . 

This implies that the Fourier coefficients satisfy the 
relation 



a 

n 



(F) 



+ nb 

n 



(F) 



0 



(15) 



The extensional mode is defined to be geometrically orthogonal 
to the flexural mode so that 



na - b 

n (E) n (E) 



= 0 



(16) 



Substituting equation 15 into equations 14 yields the result 



2 

[ pA ( (— — v jr~— ) + zn 2 ) ] a + [^z(n 2 -l) 2 ] a^ 

n 2 n (F) R 2 n (F) 



= c 



n 



) = -a /n 

n (F) n (P) 



( 17 ) 



where z = (r/R) 
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If equation 16 is substituted into equations 14, the result 
is 



[pA ( (n +1) + zn ) ] a 



n 



(E) 



r EA 

V 



2 2 2 2 
( (n + 1) + z (n - 1) ]a 



n 



(E) 



c 

n 



(18) 



b 

n 



(E) 



na 

n 



(E) 



Initial solutions were carried out using equations 
17 and 18 and then combining the results by 

a = a + a 

n n (F) n (E) 



b = b + b 

n n (F) n (E) 

This procedure gave results which satisfactorily predicted 
the ring bending moments but gave poor accuracy for axial 
force determinations. 

Re-examination of the foregoing solution process 
revealed that the mode shapes defined by equations 15 and 
16 are not orthogonal with respect to the mass or the 
stiffness matrices of the system. It was accordingly 
decided to return to equations 14 and solve them 
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simultaneously in the time integration algorithm. This 
change led to accurate axial force determinations. The 
separated equations of motion, 17 and 18, are used in the 
program as described in Section III, b. 
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III. TIME INTEGRATION METHOD 



A. CENTRAL DIFFERENCE ALGORITHM 

Time integration is accomplished using a central 

difference algorithm. If a , b , and c are known at 

n n n 

time t^, then a^ (the value of a at t^) may be 

n n 

evaluated from the first of equations 14. Letting h 

. ( i — is ) • 

represent the time step and represent a n at time 

( i ) * 

t - h/2, the next value of a n is calculated from 



•tt+H) 

ci 

n 




+ ha 



(i) 

n 



(19) 



Using this result, the next a n (at time t^ + ^ = t^ + h) 
is calculated from 

a (l+1) = a (l) + h a (l+i5) (20) 

n n n 

** ( i ) 

The value b n is similarly found from the second of 

equations 14. A pair of equations paralleling 19 and 20 

are used to calculate b^ +2 ^ and b^. When these steps 

n n 

have been completed for each n, the value of radial 
displacement w is found at each structural node from 
equation 4. These displacements are passed to the fluid 
program where they furnish required interface boundary 
conditions to allow an advance to time t^ +1 ^. Values of 
fluid pressure at the interface nodes are returned by the 
fluid program. 
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The new nodal pressures are used to calculate c n 's at 

t^ + ^. It is then again possible to advance a n and a n 

using equations 14, 19, and 20 and to perforin the parallel 

calculations for b and b . 

n n 

The solution process is started with the ring at rest 

under uniform hydrostatic pressure. Under the loading a^ 

• ( 0 ) 

is the only nonzero Fourier coefficient. Because a = 0 

n 

. (-%) 

is known (rather than a^ , a fictitious starting value of 
• (-k) 

a is first calculated from 

n 




h •• ( 0 ) 
2 a n 



A similar starting procedure is used for t> n . 

B. SELECTION OF THE NUMBER OF 
VIBRATIONAL MODES USED 

The number of vibrational modes that can be modeled 
accurately is limited by the numerical integration algorithm. 
Specifically, the algorithm becomes unstable for time steps 
in excess of about 0.3 of the period of the structural mode. 
The accuracy of the algorithm deteriorates even before the 
stability limit is reached. For this reason, a criterion 
was established for limiting the number of modes used, 
based on the time steps selected. 

In the program the separated equations of motion for 
the extensional and flexural modes were used to find the 
frequencies needed in applying the criterion. From 
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equations 17 and 18 the natural frequencies of the 
extensional and flexural modes were found to be 




2 




n = 0,1, 2, 3 



n (F) 



2 




Cx 2,r 2 n 2 (n 2 - 1) 2 
R ; 




n 



2,3,4 



, • ♦ • 



The criterion used was that the period for the highest 
mode retained be at least five times the time step used. 

Since the extensional modes have higher frequencies for a 
given n, this criterion comes into effect first for the 
extensional modes. Once the point is reached (where 
n = n„(max)), the program switches from the coupled equations 
of motion 14, to the separated equation 17. When this 
occurs the higher extensional mode coefficients (n > n„(max)) 

ill 

are not needed, and only the flexural mode coefficients 
(for n > n„ (max) ) are used. If the time step is large 

ill 

enough that the criterion is also met by the flexural 

modes (n = n_ (max) ) then both the a 's and b 's are omitted 
F n n 

for all modes where n is greater than n„ (max) . 

r 
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IV. RESULTS 



To check the accuracy of the results obtained from the 
program, two kinds of comparisons were made. First, the 
solution of equations 14 arrived at by the program were 
compared to an analytic solution of equations 14. Second, 
to verify that equations 14 correctly predict ring behavior, 
calculated dynamic response of the ring to two particular 
loadings was compared to known static results [Reference 6] 
for the same loadings . 

A. COMPARISON OF MODEL RESULTS TO EXACT SOLUTION 

The check on solution accuracy was based on the exact 
solution for a ring initially at rest and suddenly loaded 
by a steady pressure 

P = [1.625 - .25 cos 0 - .375 cos 20] MPa 

The exact solution, using the parameters specified in 
Table I, is 



a 



a 



1 



0 



4.05844x10 3 [1 - cos (1011. 303t) ] 

-159 . 483t 2 - 1 . 5617xl0~ 4 [1 - cos (1429. t) ] 



a 



2 



- . 104129 [1 - cos (85 . 616) t] 

- 3. 7355xl0~ 5 [1 - cos (2259. 7t) ] 



( 22 ) 
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Table I gives results for radial displacements at 
t = 1 msec for 8=0°, 90° and 180°. Computer results are 
given for four different time steps. Even the coarsest 
time step gives results within 1% of the exact values 
obtained from equations 22. 

B. COMPARISON OF MODEL RESULTS TO STATIC BEHAVIOR 

To compare the dynamic response predicted by the program 
to the known static response of the ring for the same 
loading, two simple cases of loading were assumed. The 
first loading case examined was uniform pressure surrounding 
the ring; the second case was two equal and opposite concen- 
trated loads acting 180° apart as shown in Figure 3. The 
second loading case could not be represented exactly by 
the finite trigonometric series loading representation used 
by the program. It was approximated by nonzero values of 
pressure only at 9 = 0° and 8 = 180°. 



B 




Figure 3. Concentrated Loading on Ring 
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TABLE I 





Time 

Step 

sec 


AR 

NP1 

m x 10 7 


AR 

NP2 

m x 10 7 


AR 

NP3 

m x 10 7 


AR 

predicted 
by program 


2.5xl0~ 4 


11619. 


23687. 


17556. 


IxlO -5 


11679 


23472. 


17553 


lxl0~ 6 


11682 


23469. 


17555. 


Lxl0~ 7 


11683. 


23469. 


17555. 


AR from 
equations 17 




11683. 


23469. 


17555. 



2 

Ring Parameters: R = 5 m, A = .05 m , E = 200 GPa 

p = 7830 kg/m^, r = .158m 



Comparison of program solution with 
analytic solution 



TABLE II 





AR (max) 
mm 


Axial 

Force 

(Max) 

MN 


Program results 


4.98 


-9.8 


Twice Known 


5.00 


-10.0 


Static Results 







2 

Ring Parameters: R=5m, A = .05 m , E = 200 GPa 

p = 7830 kg/m^, r = .158 m 



Comparison of program dynamic solution 
with twice known static solution for 
uniform external pressure on ring 
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Comparison between the dynamic and static responses 
of the ring were made for the radial displacements, bending 
moments, and forces at particular nodal points on the ring. 
For the uniform pressure case any nodal point on the ring 
can be used since the displacements are uniform over the 
ring. For the concentrated loading case, the nodal points 
chosen were at the points of application of the forces and 
midway between the application points. 

For the dynamic response of the ring predicted by the 
program, the ring is taken to be initially at rest and 
undeformed. The anticipated behavior of the ring after 
the application of the loading is for the ring to deflect 
through the radial displacements found for the static 
loading and to continue to deflect until a radial displace- 
ment with a magnitude of approximately twice the static 
displacement is reached. The maximum deflections were 
expected to occur after a time approximately equal to 
one half the period of the vibrational mode that dominates 
the deflection of the ring. For the uniform loading case 
the dominant vibrational mode is the fundamental extensional 
mode. For the concentrated loading case, it is the second 
flexural mode. The periods for these vibrational modes can 
be found using equations 21. Tables II and III give the 
comparisons between the dynamic response from the program 
and twice the known static response for the two loading 
cases considered. As is seen, relatively good correlations 
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TABLE III 



dumber of 
Stodal Points 
Used 




AR^(max) 

m 


ARg (max) 
m 


Axial 
Force 
at B 

MN 


Bending 
Msment 
at B 

MN*m 


cr 


Model 

Results 


.2965 


-.2839 


-3.760 


-7.910 


j 


Static 
Analysis 
Results 
Times 2 


.2634 


-.2421 


-3.536 


-6.434 ' 


7 


Model 

Results 


.2013 


-.1857 


-2.353 


-4.737 


Static 
Analysis 
Results 
Times 2 


.1862 


-.1713 


-2.500 


-4.543 


Q 


Model 

Results 


.1510 


-.1393 


-1.743 


-3.544 


j 


Static 
Analysis 
Results 
Times 2 


.1425 


-.1311 


-1.913 


-3.477 



2 

Ring Parameters: R = 5 m, A = .05 m , E = 200 GPa 

p = 7830 kg/fai^, r = .158 m 



Comparison of dynamic solution with twice known 
static solution for loading of Fig. 3 
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between the dynamic responses and twice the static responses 
were obtained. The maximum values shown in Table III occur 
at approximately 35.3 milliseconds after time zero, which 
compares favorably to one half of the period of the second 
flexural mode (36.6 milliseconds). The maximum values in 
Table II occur at 3 milliseconds which is approximately 
one half the period of the fundamental extensional mode 
(3.11 milliseconds). 

C . OUTPUT EXAMPLES 

A sample of the output produced by the structural 
program is shown in Figures 4 and 5. The print code at the 
top of Figure 4 is an input code used to select the informa- 
tion desired in the print out. Following the print code 
are two sets of input parameters to be used in the structural 
model. These are followed by the mass and stiffness 
coefficients calculated for each of the vibrational modes. 
This information remains constant for any given run and 
is therefore only printed once. The rest of the information 
varies with time and can be printed as often as desired. 

The pressures at nodal points generated by the fluid model 
are printed out; these pressures are for one time step 
before the time printed below the pressures. Following the 
time come the Fourier coefficients for the radial displace- 
ments, tangential displacements, and pressures. In 
Figure 5 at the top, the bending moment and the axial force 
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at each nodal point are printed for the time given in 
Figure 4. Tabulation of the radial displacement follows. 

Figure 6 is an example of the graphical output obtained 
from the fluid model when the fluid and structural models 
are run jointly. It is a time sequence of eight plots of 
fluid nodal pressures over the domain. The plots are shown 
at 8 ms intervals. The left hand edge of each plot repre- 
sents the plane of symmetry and the top and bottom rows 
represent, respectively, the entry and exit faces for the 
shock. The x's on the lower left side of the plots 
represent dummy nodes inside the structure. The pressure 
ranges for the mapping characters are shown in Figure 7. 

In Figure 6 the development of the cavity can be followed. 

It develops in the second, third, and fourth frames, 
collapses in the fifth, and has vanished by the sixth. 



34 



PRINT CODE : 121111111 1 112 1 141512 1712 1 1 
NNEL NNP M N NNPR1N NNPP2N 
2 @ '? 17 



16 17 



R 



0 . 50OD-O 1 0 . 200 B 12 



. yy 



o 

9 

10 

11 

12 

13 

14 

15 

16 
17 



XKfllN 

0.4004OD 09 
0.40000B 09 
0. 40366D 09 
0. 42560D 09 
0. 4900OD 09 
0.63040B 09 
0. 390O0D 09 
0. 13216D 10 
0. 19376D 10 
0.29600D 10 
0. 43 20 4 D 10 
0.616O0D 10 
0. 25796 D 10 
0. 11690D 11 
0. 15610B 11 
0.20470D 11 
0. 26410D 11 



KKR2N 

0 . 0 

0.10217H 07 
0. 20424-0 07 
0.3065 1C 07 
0.40S6SD 07 



0.5102 



0.61303D 07 



0.715 

0.2172 



6D 



0B 07 
7D 07 



0.91954D 07 
0. 10217D 03 
0.11239H 03 
y. 1226 ID 02 
0. 1 322213 03 
0. 14204D 03 
0. 15326B 08 
0. 16347D OS 



PI 

0. 158 

KMR1N 
391 . 50 
391 . 89 
393.07 
395. 02 
397. 76 
401.29 

405. 59 
410.63 
416.56 
423.21 

430.65 
433. 37 
447. 33 

457. 66 
463. 23 

479.59 
491.72 



MODAL POINT PRESSURES 



0. 12D 08 


0. 98 D 


07 0.56B 07 


0.21B 07 i 


0.2 ID 06 


0.20D 


06 0 . 20 D 06 


0.20B 06 1 


0. 16D 06 


0.20D 


06 0. 20 D 06 


O.20D 06 1 


TIME = 


0. 250 


MS 




BN 




AM 


CM 


0. 0 




0. 599 23 B- 03 


0. 162973 


— M . *i!j't.*--6'D 


-05 


0. 19121D-03 


y • ui' -5 »■ jl 


— 0. 76030D 


-05 


0. 16725B-03 


0 B 237373 


-0.90952E 


—65 


L ! • 1 o4c5 j.--— 03 


0 ! 1839 ID 


-0. 8786 2D 


-Q5 


L-1 m *— 4 J— f 1 — ^' 1 4 


0. 12588D 


-0. 70538D 


-05 


0.6391 SB-04 


C. 86^200 


y • 4868 1 D 


-05 


0.37398D-0- 


y. 4!~ yy cD 


-y. 2781 2D 


-05 


0. 187321-0- 


0.221721! 


-0. 1 1520D 


-05 


0.6959 ID-05 


76529. 


-0. 2-765D 


-y6 


0 . 22223D-05 


1S334. 


0. 7676 2D 


—03 


—0 .767 62E-07 


-8451.9 


0. 705 17D 


-07 


-0.77569D-06 


-15526. 


0. U450D 


-07 


-0. 137-OD— 06 


-4325. 4 


0. 64630D 


-07 


-0 . 340 1 SB-06 


— 1 4 c 4 c • 


0.4838 ID 


-07 


— y . 68488'D— y6 


-10763. 


8. 38325B 


-07 


-0 . 57 4870-06 


-8137. 1 


0.3383 ID 


-07 


-0. 541000-06 


-9139. 7 



RHO 

0..733B 04 
XKE1N 

0 . 0 

0.40000D 09 
0 . 30 000 D 09 
0. 1 280-0 D 10 
0. 16000D 10 
0. 20000 D 10 
0.24000D 10 
0.23000D 10 
0. 32Q00D 10 
0. 36000 D 10 
0. 40000 D 10 
0 . 44000. D 10 
0. 48000 D 10 
0.520001) 10 
0. 56000 B 10 
0. 60000 D 10 
0.64000D 10 



0.52B 06 
0.21 D 06 
0. 20 B 06 



H . c’JD 06* 
0.2 IB 06 



Fig. 4. Output Example 
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IMP MOMENT FORCE 





1 


6 1 *“?*bb • 


-S.33717D 




C 




-0.31430D 




* 


-31*43. 






4 


-46163. 


- 0 ! 1 4334 B 




cr 

.j 


-37598. 


—0 . 1 1 448 D 




6 


-7468. 4 


-0. 18173D 




—7 

i* 


-5081 . 3 


-0 . 99940 B 




i"i 


-FI005. 4 


-0. 10845D 




Q 


-5153.0 


-0. 18006D 




10 


-4923. 6 


-0 . 998S5B 




11 


-5337.5 


-0. 10017D 




13 


-3399. 1 


-0. 10010D 




13 


-7444.3 


—0 . 997931- 




14 


-5635. £ 


-0.99765D 




15 


-5365. 8 


— 0 i. 18009D 




16 


—4963. 8 


-8. 99989 D 




17 


-4983.9 


-0.99886D 


HP 


NODAL PO 


I NT DISPLACEMENTS 




1 




0 . 13157D-03 




•-» 

U. 




0 . 1 1784D-02 




O 




0 . 86677D-03 




4 




0.63037D-03 




5 




0. 51S33D-03 




ti 




0 . 50099 D— S3 




i’' 




0 . 50035B-03 




3 




0 . 49977D-03 




9 




0. 49981D-03 




10 




0. 49984D-03 




11 




0 . 50833D-03 




1 £ 




0. 50019D-03 




13 




0 . 49695D-03 




14 




0 . 49939D-03 




15 




0. 49947D-03 




16 




0. 49931D-03 




17 




0 . 49944 D— 03 




T= 


10. 40-- 13. 7 


3 l0.c8.53 





07 

87 

y i*' 

87 

87 

87 

8 b 

87 

87 

0 b 

07 

87 

8b 

0b 

87 

0b 

06 



Fig, 5. Output Example (cent.) 
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CCCCCEEEEEEmE££EBH£E£EECCC£ECB 
CCCCCE EICCCCCCCCCCCCCCCCCCCCCCCCC 
E’DCCCCCC'! I'DBrLrDtiriDDr'BDPr'rDDDDCCC 
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T - 36.000 « 
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Symbol 



Pressure Range (m pa) 



X 

blank 

1 

3 

6 

8 

H 

A 

B 

C 

D 

E 

F 



Structure 
Cavitated 
0 - .04 
.04 - .09 
.09 - .14 
.14 - .18 
.18 - .43 
.43 - .77 
.77 - 1.10 
1.10 - 1.43 
1.43 - 1.77 
1.77 - 2.10 
over 2.10 



Figure 7. Pressure Code 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The comparisons made between the program results where 
convergence of the solution has been obtained to results 
obtained by other methods, show that the model satisfactorily 
predicts the dynamic response of a ring to external pressures. 
The ring model cannot predict exactly the behavior of a real 
structure such as a submarine hull. However, by a suitable 
choice of ring parameters the response of the ring structure 
will provide a good first approximation to the response of 
a real structure. 

A useful improvement to the program would be to include 
a section which would automatically pick out maximum values 
of bending moments and axial forces in the ring as well as 
when and where they occur. This would allow any desired 
time intervals between data printout to be used without 
missing peak responses. 



39 



APPENDIX 



PROGRAM LISTING 
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